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We study possible factorizations of supersymmetric (SUSY) transformations in the one- dimensional quan- 
tum mechanics into chains of elementary Darboux transformations with nonsingular coefficients. A clas- 
sification of irreducible (almost) isospectral transformations and of related SUSY algebras is presented. 
The detailed analysis of SUSY algebras and isospectral operators is performed for the third-order case. 

1. Introduction: definitions and notation of the SUSY QM 

The concept of supersymmetric Quantum Mechanics (SUSY QM) represents an algebraic form of transfor- 
mations of (complete or partial) spectral equivalence between different dynamical systems [I]-[S]. At present, 
there is a number of reviews [5]-[T2] devoted to development and various applications of the SUSY QM; the 
reader is referred to these reviews, which are addressed for a more detailed study of this approach to construction 
of isospectral systems. Isospectral transformations of that kind are the Darboux-Moutard-Crum transforma- 
tions [13j . |14|-[T7]. which are are known in the theory of ordinary differential equations for a long time 0- 
In the simplest cases, intertwining of two differential operators (for instance, Hamiltonians of one-dimensional 
quantum systems) by means of Darboux operators entails their factorization into differential multipliers which 
are formed by the same Darboux operators (Schrodinger factorization [18j, [19] and its generalizations [20] . 
[21]). However, in general, this is not the case, and both interrelation between pairs of dynamical operators 
("Hamiltonians") with (almost 0) equivalent spectra and structure of operators which generate the spectral 
equivalence are not that simple [22] - [IS] . Precisely this interrelation in the one-dimensional QM is the focus of 
the present paper. In particular, we present rigorously justified answers to the following questions: in what 
cases can the higher-order Darboux-Crum transformations be constructed with the help of a sequence of inter- 
twining transformations of lower order which relate a chain of (almost) isospectral intermediate Hamiltonians 
with real nonsingular [f| potentials; what are elementary blocks for a nonsingular factorization of intertwining 
operators ; in what way is the irreducibility of elementary blocks of isospectral transformations indicated in 
the SUSY algebra and in the structure of kernels of those transformations? The structure of the paper is as 
follows. After a short reminder of notation and basic definitions of SUSY theory of isospectral transformations 
we formulate basic theorems on the structure of a polynomial SUSY algebra and on minimization of this al- 
gebra up to its essential part (proofs of these theorems can be found in our preceding paper [26]). Then we 
present a classification of irreducible (almost) isospectral transformations and related SUSY algebras (partially 
described in [T^j, [17]-[3T]). Next, we define a potential class K that is invariant under transformations of the 
Darboux-Crum type and formulate two theorems on reducibility of differential operators of spectral equivalence 
transformations. The paper is completed with a detailed analysis of the third-order SUSY algebras and isospec- 
tral operators as a first stage in proving above-mentioned theorems on reducibility. A complete proof will be 
published in a forthcoming issue. 

Let us start with a definition of the SUSY algebra and notation of its components. Consider two one- 
dimensional Hamiltonians of the Schrodinger type — —d^ -\- Vi{x) and /i^ = -f V2{x), d = d/dx, which 

Translated from Zapiski Nauchnyikh Seminarov POMI, Vol.335, 2006, pp. 22-49 . 
^In the monograph 1171 the Darboux transformations are given for a wider class of partial differential equations including 
non-stationary Schrodinger one and some nonlinear equations. 

■^We say that operators have almost equivalent spectra if their spectra are different only at a finite number of eigenvalues. 
^In this case, the potentials are sufficiently smooth, but potentials having singularities weaker than 1/x^ are also acceptable. 
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are defined on the entire axis and have nonsingular potentials Vi,2{x) ■ We assemble the Hamiltonians into a 
super-Hamiltonian, 

(";,.")■ ™ 

Assume that the Hamiltonians and h~ have an (almost) equal energy spectrum of bound states and equal 
spectral densities of the continuous spectrum part; let such an equivalence be provided by the Darboux-Crum 
[131 114j operators with the help of intertwining, 

h'^QN = 1n>^^ = h^lN- (2) 

Further on, we restrict ourselves to differential Darboux-Crum operators of finite order N , 

N 

4-E^fc(^)^'' ^%^{Tlf, (3) 

fe=0 

with real, sufficiently smooth coefficients w^{x). In this case, in the fermion number representation, the non- 
linear M — \ SUSY QM is formed by means of nilpotent supercharges, 

Qa^ = ( f ) , ^^={q-^ ) ' Qn = QI = 0. (4) 
Obviously, the intertwining relations ^ lead to the supersymmetry of the Hamiltonian H , 

[H, Qn] - [H, Qn] - 0. (5) 
This nonlinear SUSY algebra is closed by the following relation between the supercharges and Hamiltonian, 

{Qn,Qn} = Pn{H), (6) 

where Piq{H) is a differential operator of 2A^th order commuting with the Hamiltonian. Depending on a 
relation between the supercharges Qn,Qn (the intertwining operators g^), the operator Pn{H) can be either 
a polynomial of the Hamiltonian if the intertwining operators are connected by the operation of transposition: 
— [qj]) = X]feLo(~^)'^'^fe (*^)' '-'^ ^ function of both the Hamiltonian and a differential symmetry operator 
of odd order in derivatives (see a detailed analysis and references in |26)). In our paper, we confine ourselves 
with the first case in which the conjugated supercharge is produced by transposition, Qn = Q% (a relevant 
theorem on the structure of such a SUSY is formulated below). 

2. Basic theorems on the structure of QM with a nonUnear SUSY 

Theorem 1 (on supersymmetric algebra with transposition symmetry). 

Let (f)f^(x), n — 1, . . . , N be a basis in keiq^: 

cL^4>t = 0, = (g+)*. (7) 

Then: 

1 ) the action of the Hamiltonians on the functions (fx^ [x) is described by constant N x N matrices, 

N 

h^4>t^Y.^nm4>^ n = l,...,iV; (8) 

rn — 1 

2) the closure of the supersymmetry algebra takes a polynomial form, 

{Q, Q*} = Aet[El - S+]e=h = det[£;i - S-]e=h = Pn{H), (9) 

where I is an identity matrix and is the matrix with entries S^^. 
Corollary 1. The spectra of the matrices S+ and S~ are equal. 
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In what follows, for an intertwining operator, its matrix S is defined as the matrix which is related to 
operator in the same way as are related to . In this case, we do not specify the basis in the kernel of 
the intertwining operator in which the matrix S is chosen if we concern only with spectral characteristics of the 
matrix or, that is the same, spectral characteristics of the restriction of the corresponding Hamiltonian to the 
kernel of the intertwining operator considered (cf. (H))). 

A basis in the kernel of the intertwining operator in which the matrix S of this operator has a Jordan form 
is called canonicaZ; elements of a canonical basis are called transformation functions. 

Assume that the intertwining operators are represented as a product of the intertwining operators k^_j^,j 
and p^, < Af < so that 

9n = Pm^n-m^ 9n — ^n-mPm'i Pm^m = f^^Pti- Pm^^ = ^mPm'^ 

^N-M^^ = hMk'lj_j^j, k'^^M^M = h^k'^_i^_j, hM = -d^ + vm{x), (10) 

where the coefficients k^_j^,j and p^j as well as the potential vm{x) may be complex and/or singular. The 
Hamiltonian Hm is called intermediate with respect to /i"*" and . In this case, by Theorem 1, the spectrum 
of the matrix S of the operator is a union of the spectra of the matrices S for the operators k^_j^j and pf.j. 
The potentials Vi {x) and V2 {x) of the Hamiltonians and h~ are interrelated by the equation 

V2{x)^Vi{x)-2[lnW{x)y', (11) 

where W{x) is the Wronskian of elements of an arbitrary (a canonical as well) basis in kerg^. The validity of 
Eq. pT|) follows from the Liouville-Ostrogradsky relation and the equality of coefficients at in q]^h^ and 
h~qjj (see the intertwining in 

An intertwining operator q^ is called minimizable if this operator can be presented in the form 

g± =P(/i±)p± =p±P(/.T), (12) 

where p^j is an operator of order M which intertwines the same Hamiltonians as ( i.e. p^jh^ — h^p^j) 
and P{h^) is a polynomial of degree (A^ — M)/2 > 0. Otherwise the intertwining operator q^ is named as 
non-minimizahle. 

The following theorem contains necessary and sufficient conditions under which an intertwining operator is 
minimizable or not (a proof can be found in [26] ). 

Theorem 2 (on minimization of an intertwining operator) 

An intertwining operator can be presented in the form 

m 

QN=PMl[^^i-h^)""^ (13) 

1=1 

where p^j is a nonminimizable operator intertwining the same Hamiltonians as q^ (so that p^jh^ = h^p^j), if 
and only if a Jordan form of the matrix S of the operator q^ has m pairs ( and no more ) of Jordan cells with 
equal eigenvalues Xi such that, for the l-th pair, Ski is the order of the smallest cell and ki + dki is the order of 
the largest cell . In this case, M — N — 2J2]^i = X]r=i where the ki,m + l^l'^n are orders of the 
remaining unpaired Jordan cells. 

Remark 1. A Jordan form of the matrix S of the intertwining operator q^ cannot have more than two 
cells with the same eigenvalue A; otherwise ker(A — h^) includes more than two linearly independent elements. 

Corollary 2. Jordan forms of the matrices S of the operators q^ and q'^ coincide up to permutation of 
Jordan cells. 

If a Jordan form of the matrix S of an intertwining operator has cells of order higher than one, then the 
corresponding canonical bases contains not only formal solutions of the Schrodinger equation but also formal 
associated functions, which arc defined as follows [321. 

A function ^n.i{x) is called a formal associated function of i-th order of the Hamiltonian h for a spectral 
value A„ if 

(/i - A„)^+Vn,» = 0, and (/i - A„)Vn„ ^ 0. (14) 

The term 'formal' emphasizes that this function is not necessarily normalizable (not necessarily belongs to 
i2(R))- In particular, an associated function tl^n.o of zero order is a formal eigenfunction of h (not necessarily 
a normalizable solution of the homogeneous Schrodinger equation). 
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3. Classification of really (ir)reducible SUSY transformations 



The intertwining operator q^^ is called {really) reducible if this operator can be presented as a product of two 
nonsingular intertwining operators (with real coefficients) k^_M ^^idp^, < M < so that Eqs. pU]) are valid 
and the intermediate Hamiltonian Km has a real nonsingular potential. Otherwise is called {really) irreducible. 

Really irreducible, nonminimizable, intertwining operators of second order with real coefhcients can be 
divided into three types [12] . 

A really irreducible intertwining operator of I type is a differential intertwining operator with real coefhcients 
for which eigenvalues of the matrix S have nontrivial imaginary parts and are mutually complex conjugate. 

Let us show that any intertwining operator satisfying this definition is, in fact, really irreducible (the 
case of q2 is treated similarly). Indeed, let Lp'^2{^) be a canonical basis of kerg^ such that — M,2'Pi2i 

\l ~ X2 ^ Xi. Assume that q2 is reducible, i.e., there exist intertwining operators and p^ with real 
nonsingular coefficients such that 

Q2=KPi^ Pih^ = hiPi, k^hi^h^k^, (15) 

where hi is an intermediate Hamiltonian with a real nonsingular potential. Obviously, a basis in the kernel of 
Pi consists either of ip^ or of ip2 ■ We restrict ourselves to the case of since the case of ip2 can be considered 
in the same manner. Then p^ = d — {ip^)' /ip^ , and, consequently, f{x) = {pf^)' jpf^ is a real- valued function. 
But then 

Pl (x) = Ce^ C = Const; 

hence, 

^i(^)-Ai = ^^^ = /^(x) + /'(.) 

{X) 

is a real- valued function as well, and we get a contradiction with the condition that 7^ Ai. Thus, any operator 
that satisfies the above definition is indeed really irreducible. 

The degenerate case V2,\{x) = Const should be singled out. In this case, = = hi, and the canonical 
basis ker q^ can be chosen in the form 

p>t{x)=e'^\ ^±(a;)=e^*^ k ^ k* (16) 

so that eigenvalues of the matrix S of the operator qf and the operator itself are as follows: 

h^ipf 2^ >^l2ft,2^ AjF=F24-fc^ AJ=F24-fc*^ qf ^d^ -2Rekd+\k\^. (17) 

Note that potentials of the intermediate Hamiltonians which correspond to two possible factorizations of a 
really irreducible intertwining operator of the I type into intertwining operators of first order, i.e., 

F2,i(x)-2[ln^±2(x)]", (18) 

where '/'^2(^) ('Pr2(^)) ^ canonical basis in kerq^ (kerq^), always have a nontrivial imaginary part (see [23j l 
with the only exception of the case V2,i(a;) = Const. 

A really irreducible intertwining operator of the II type is a differential intertwining operator q^ of second 
order with real coefficients such that: 

(1) eigenvalues of the matrix S of the operator qf are real and different; 

(2) both elements pf{x) and p^{x) of a canonical basis of kerg^ have zeroes. 

The irreducibility of intertwining operators satisfying this definition follows from the fact that otherwise the 
equalities, 

q+=pfk+, p+hi^h+p^, k+h-^hik+ (19) 

take place, or, according to a basis in kerfcj*' (kerp^) consists either of pf{x) or of p}f{x), and the 

potential of the intermediate Hamiltonian hi is described by one of Eqs. psp . i.e., has a singularity(ies) by the 
second item of the definition. We also note that potentials of intermediate Hamiltonians which correspond to 
two possible singular factorizations of a really irreducible intertwining operator of the II type into intertwining 
operators of first order given by (|18p are real since the pf 2i^)'-^^^ be always chosen real. 
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A really irreducible intertwining operator of the III type is a differential intertwining operator q2 of second 
order with real coefficients such that: 

(1) the eigenvalues A1.2 of the matrix S of the operator are equal, Ai = A2; 

(2) a canonical basis in kerq^ consists of formal eigenfunctions, ipfQ{x), and associated functions, (pfi{x), of 
the Hamiltonian which assemble into a Jordan cell, 

(3) <pfo(^) least one root. 

The irreducibility of an intertwining operator satisfying this definition follows from the fact that otherwise 
equalities ([TO)) take place, or, according to (fT5|) . a basis in kerfcj*" (kerp^) consists of ipfQ^x), and a potential of 
the intermediate Hamiltonian hi is described by the equation 

V2Ax)-2[ln^%{^)]", (20) 

i.e., has a singularity (ies) by the third item of the definition. The potential of the intermediate Hamiltonian, 
which corresponds to the only possible singular factorization of a really irreducible intertwining operator of the 
HI type into intertwining operators of first order given by ([^0)) . is real since the ififoix) can always be chosen 
real. 

Obviously, other types of really irreducible nonminimizable intertwining operators of second order do not 
exist. 

Further on, we formulate two assertions which characterize reducibility of intertwining operators of any order 
in an exhaustive way: 

assertion (1) of Theorem 3 on the reducibility of a nonminimizable intertwining operator with real spectrum 
of the matrix S, multiplied by an appropriate polynomial of the Hamiltonian, into (a product of) intertwining 
operators of first order; 

assertion (2) of Theorem 4 on the reducibility of a nonminimizable intertwining operator with arbitrary 
spectrum of the matrix S into (a product of) intertwining operators of first order and irreducible second-order 
intertwining operators of the I, II and HI type. 



4. Theorems on complete reducibility of intertwining operators 

In what follows, we use a class K of potentials V{x) such that: 

1) V{x) is a real- valued function from C^; 

2) there exist numbers i?o > and e > {Rq and e depend on V{x)) such that the inequality V{x) ^ e 
takes place for any > Rq; 

3) the functions 



±Ro 



are bounded for x ^ Rq and x ^ —Rq, respectively. 

In addition, we discuss normalizability and nonnormalizability of functions at -|~oo and/or at —00; these 
properties are defined as follows. 

A function f{x) is called normalizahle at +00 (at —00) if there exists a real number a+ (a_) such that 



+00 

2 



f{x)\^dx<+oo / |/(x)|^(ix<+oo . (22) 



Otherwise /(x) is called nonnormalizable at +00 (at —00). 

Theorem 3. (on reducibility of "dressed" nonminimizable intertwining operators) 

Assume that the following conditions are satisfied: 
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1) h+ = —d^ + Vi{x), Vi{x) G K, and the potential V2{x) of the Hamiltonian h~ is real and continuous; 

2) h+ and h^ are intertwined by a nonminimizable differential operator of Nth order with coefficients 
from Cg, so that 

q^h+ = h^q]^; (23) 

3) the algebraic multiplicity of Xi, the ith eigenvalue of the matrix S for the operator q~^, is equal to ki, 
i = 1, . . . , n, so that fei + • ■ • + fc^ = iV; all of the numbers Xi are real and satisfy the inequalities 

^ Ai > A2 > . . . > A„; (24) 

4-) is the spectrum of the matrix S of the operator q^; 

5) Ei±, i — Q, 1, 2, . . . , is the energy of the ith (from below) bound state of h^ ; is the number of bound 
states of /i^ with energies of which are included into A; N± be a number of bound states of h^ with energies 
not exceeding Ai; 

6) 

P±{E)^ W {E-E,±). (25) 

Ei±<\i,Ei±(^A 

Then: 1) V2{x) G K; coefficients of q~^ belong to and are real; = ('Zjv)* ^'^^ '''^^^ coefficients from 
and intertwines h+ and h~ , so that 

h+q+^q+h-; (26) 

2) P+{E) EE P-{E); the degree of P±{E) is equal to N+ - N+ = - N- ; 

3) the operator q^P±{h^) intertwines h+ andh~ and can be presented as a product of N+N++N^—N+—N~ 
intertwining operators of first order with real coefficients from , so that: 

a) potentials of all the intermediate Hamiltonians belong to K ; 

b) the eigenvalue of the matrix S of the l-th operator (from the right) in the factorization under consideration 
is equal to Ei-i^±, 1 = 1, . . . , N± and an element of the kernel of this operator is normalizable at both infinities; 

c) the eigenvalue of the matrix S of the l-th operator (from the left) in the factorization under consideration 
is equal to Ei^i^z^, 1 = 1, and an element of the kernel of this operator is nonnormalizable at both 
infinities; 

d) the set of eigenvalues of the matrices S for operators from the N± + l-th to the N± + N — N+ — -th 
one (from the right) in the factorization under consideration coincides witl^ A \ {{Em^} U {i?i_}). In addition, 
the eigenvalue of the matrix S for an operator of this group does not decrease as the number of the operator 
increases (from the right to left); a basis element of the kernel of any operator in this group is normalizable at 
one of the infinities only. 

Theorem 4. (on complete reducibility of nonminimizable intertwining operators) 
Assume that the following conditions are satisfied: 

1) h+ = —d^ + Vi{x), Vi{x) € K; the potential V2{x) of the Hamiltonian h~ is real and continuous; 

2) h+ and h^ are intertwined by a nonminimizable differential operator qj^ of Nth order with real coefficients 
from C^, so that 

q]^h+ = h^q'^] (27) 

3) the algebraic multiplicity of Xi, the ith eigenvalue of matrix S for operator qj^, is equal ki, i = 1, . . . , n, 
so that fci + • • • + fc„ = N ; the set of values Xi contains M real values and L pairs of mutually complex conjugate 
ones, so that M + 2L = n; the numbers i = 1, . . . , M correspond to real Xi, and A,; > Ai+i, i = 1, ■ ■ ■ , M — 1; 

4) if Xi is real, then Ai ^ 0; 

5) Ei±, i = 0, 1, 2, . . . is the energy of the ith bound state (from below) of h^ ; K± = ma.x{i : A,; > Eo±}, if 
Ai > Eq±, and K± = 0, if either Ai ^ Eq± or ImAi ^ 0. 

Then: 1) V2{x) G K; coefficients of qjj belong to ; = {q)^Y has real coefficients from and 
intertwines h+ and h^ , so that 

h+q+=q+h-; (28) 

^In this formula, one has to take into account multiplicities of eigenvalues as follows: if A is contained in A with algebraic 
multiplicity Ki, in {Ei^} with multiplicity K2 and in with multiplicity Ki (obviously, K2 and K3 can take values and 1 

only), then the value A is contained in A \ {{Ei^} U with multiplicity Ki — K2 — K3 if Ki > K2 + K:j or is not contained 

if <iK2-\-K3. 
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2) can be presented as a product of really irreducible intertwining operators of first and second order with 
real coefficients from , so that: 

a) potentials of all the intermediate Hamiltonians belong to K ; 

b ) the first 

M+L 

Ji = X! (29) 

i=M+l 

operators from the right in the factorization of under consideration have an order 2 and are really irreducible 
operators of the I type ; in addition, one can realize that the related pairs of mutually complex conjugated 
eigenvalues of the matrix S operator are ordered arbitrarily; 

c) the second (from the right) group of operators in the factorization under consideration consists of 

J2T=N- 2Ji - 2 Ja^, (30) 

operators of first order, where 

rl 1 

'^3T=[2^H' ^^^^ 

1=1 

and 

(i) if ^ ki is even, then the eigenvalue of the matrix S for the operator g^y which corresponds to the Ith 

i=l 

(from the right) of these operators does not exceed the eigenvalue related to the I + 1-th operator, I ^ \, 

(i'i') if ^ ki is odd, then the eigenvalue of the matrix S for the operator q^^ which corresponds to the Ith 

i=l 

(from the right) of these operators does not exceed the eigenvalue related to the {I + l)th operator, 1^1, . . . , 
— 2; Xk^ is an eigenvalue of the (J2=p — l)th operator and Xk^+i is an eigenvalue of the {J2^)th operator; 
in this case, the latter eigenvalue is equal to Eq^ ; 

d) the third (from the right) and the last group of operators in the factorization under consideration consists 
of J^zfi really irreducible operators of II and III type, wherein the largest of eigenvalues of the matrix S for 
the operator q^ which corresponds to the Ith of these operators (from the right) does not exceed the smallest 
eigenvalue of the matrix S for the operator qj^ which corresponds to the [I + l)th of these operators, I — 1, ■ ■ ■ , 

Remark 1. If i?o=p is not an eigenvalue of the matrix S for the operator g^^, then ^ ki is even since 

i=l 

otherwise the eigenvalue of q^q^ = Y[i=ii^^ ^ •^i)'°' the ground state wave function of h^ is negative. 
Proofs of these theorems will be published in a forthcoming issue. 

The next aim of this paper is to show rigorously that any nonminimizable intertwining operator of third 
order with real coefficients is really reducible Such a proof is a necessary stage in the study of reducibility of 
intertwining operators of arbitrary order. For this purpose, we first derive differential equations for Wronskians 
of subsets of a canonical basis of the kernel of the intertwining operator q^ of an arbitrary order A'^. These 
equations form a base of the proof of reducibility for an arbitrary intertwining operator of third order (Theorem 
5) and also can be used, for instance, to examine reducibility of intertwining operator in the general case where 
a Jordan form of its matrix S is a single Jordan cell. 

5. Derivation of system of equations for partial Wronskians 

Let 4>j{x), j — I, . . . , N he a canonical basis in kerq^ and let Aj, j = 1, . . . , be an eigenvalue of the matrix 
S for the operator q^ corresponding to the Jordan cell to which 4>j{x) is related. It is shown in [26], Lemma 1, 
that: 

^For preliminary consideration of this issue, see in 1331 134| 
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1) the intertwining operator can be presented as follows 

In ~ ^'i ■ ■ ■ ^AT' 

where the Darboux operators 



can be chosen to satisfy the equalities 



2) the following relations take place 



j = l,...,iV-l, 



3) the intermediate Hamiltonians hj, j = 1, . . . , N — 1 have the Schrodinger form: 

= -d^ + Vj{x), Vj{x) = x^jix) - x'jix) + Aj = x^j+iix) + x'j+iix) + Aj+i, 

Vi{x) = vn{x) xlii^) - x'n{x) + Atv, V2{x) = vo{x) = Xi{x) + x'i{x) + Ai, 

but in general with complex and/or singular potentials; 

4) the intertwining relations 



are realized. 

Let us introduce generalized Crum determinants 



/ = 0, ...,iV-l. 



W,{x) 



cj)N{x) (f>'Nix) 



By p4p . the following expressions for intertwining operators rj . . .r^ are valid: 



' N 



W,{x) 



</>,(x) </>;.(x) 
1 d 



therefore, 



QN-j+l 



] 



j = l,...,7V 



^ ---^ W,ix) 
Furthermore, the intermediate superpotentials Xji^) are as follows: 

^^^^^ ^ mx)/Wj+,ix)]' _ ^ W^i+iW 



.N. 



j = l,...,7V, W'^+i(x) = l. 



T4^,(x)/PF,+i(x) W,{x) W,+iix) 
To obtain differential equations satisfied by the Wronskians Wj (x) , we convert the expression 

,tW,-iix) 



1, il, ) 



W,{x) ' 



(32) 
(33) 

(34) 



(35) 

(36) 
(37) 

(38) 



(39) 



(40) 



(41) 



(42) 



(43) 



in two different ways. On the one hand, we take into account (|¥T|) and intertwinings ([55]) to show that 

m^iix) 



Aj)gj ■ • ■ gjv</'j-i • ■ ■ iNi^T-^ - Aj)0j-i = (Aj_i - Aj) ^ 



(44) 
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Obviously, equalities (I44p are valid not only if is a formal eigenfunction of but also if is a formal 
associated function of On the other hand, Eq. (|43|) can be transformed as follows. 



d- 



J+1 



( ^J-l 



\ Wj J \Wj+iJ 



where we use (|35]). Finally we obtain the equations 



For further purposes, it is convenient to introduce the functions Wj — Wj/Wj, 

w'j - Wj^2 + - "^^+2 - 2wj+i {wj ~ Wj+2) + Aj - Xj+i =0, j = 1 
in system (|46p . In addition, supplementing system (j47|) with the equation 



J = 2,...,iV. 



,iV-l, 



(45) 



(46) 



(47) 



(48) 



(i.e., the Schrodinger equation for Wn rewritten for wn) and summing up the last iV — ?i + 1 equations of the 
new system, we get the relations 



w'n + li'^+l + {Wn - Wn+l)^ + A„ - = 0, 



n = 1, 



(49) 



6. Parametric formulas for partial Wronskians 

If iV = 3 system (H71) . (|^5)) takes the following form 

w[ - 7«3 +wl - wl - 2^2(^1 - uls) + Ai - A2 = 0, 

^2+11^1- 2W3W2 + A2 - A3 = 0, 



u)'^ + wl — Vi + X3 = 0. 



Let us introduce the function 



G{x) = -[w[{x) + wfix) - Vi{x) +X1 + X2 + A3]. 



Equalities (|5D|) and (|52p imply the identity 



W2{wi - W3) = G - A2 



(50) 

(51) 
(52) 

(53) 
(54) 



In order to derive a formula which expresses in terms of G, let us compare two expressions for w'2- the 
expression, obtained by differentiation of the equality 



W2 



G- A2 
W1-W3' 



(55) 



which follows from (|54p . and the expression deduced from (|5ip after a substitution of ([55)1 into (|5ip instead of 
W2- By solving the appearing quadratic equation for w^, we come to the equality 



w,^w^ + ^' + P,(X) ^ (A - Ai)(A - A2)(A - A3) 

Z[U — A3j 



(56) 
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for a certain branch of {Cy + iPsiG) . It follows from ^ and ^ that 



= 2{G^) ■ ^ ^ 

Note that formulas (|55p. and ([57]) are valid if G{x) is different from an identical constant that equal one of 
the numbers Xj. Below, we show that the latter condition always takes place for an intertwining operator 
that cannot be stripped-off. 

Now we derive formulas which express wi in terms of G. We substitute (|56p into (|50p instead of W3 and 
obtain wi from the resulting expression to deduce that if 



then the equality 

holds; if for some interval 
on some interval, then 



[G'{x)Y + AP3{Gix))^0 (58) 



= + (59) 

2V(G')2+4P3(G) 

[G'(x)]2 + iPsiGix)) = 0, G'(x) ^ 0, (60) 



wi{x)=0 (61) 
on this interval, and the potentials Vi{x) and V2{x) are identical by ([TT|) . 



7. Smoothness of potentials and coefficients of intertwining opera- 
tors 

The following lemma indicates how smooth are the V2{x) and coefficients of the intertwining operators for 
a given smoothness of Vi (x) . 

Lemma 1. Assume that: 1) ^ —d^ + Vi 2{x), Vi{x) e C^, 3, and V2{x) e Cr; 

2) 

q^ ^d' + a{x)d^ + p{x)d + 7(2;); a{x),P{x), 7(2;) G C^; (62) 
3) q^ intertwines and , so that 

q^h+^h^q^. (63) 

Then: 

1) a{x) e C'i+\ P{x) e C^, 7(^2 G Cl-\ and V2{x) e Q; 

i/ie operators q^ and q^ = {q^ )* can he presented in the form 

q^ = + g2{x)d^ + [9'2{x) T 2.91 (x)]5 + [go{x) T g'i(x)], (64) 

where 

g2{x) = a{x), gi{x) = [a' (x) - /3(x)]/2, 50(2;) = 7(2^) + [""(2:) - P'{x)]/2; 

in addition g2{x) G Cr+\ gi{x) G Cr+\ go (a;) G C^^^ ; 
3) q^ intertwines h'^ and h~ , so that 

h+q+^q+h-. (65) 



Proof. Let us check first that 

a{x) G Cl-\ (i{x) G 7(2;) G C^-K (66) 

Indeed, inclusions ([66]) follow from relations p2p and (|40|) for j = 1 and = 3, from the fact that 4>i{x), 02 (2;) 
and (/>3(2:) belong to C^'*'^ as formal eigenfunctions and (possibly) associated functions of and from the fact 
that Wi{x), a Wronskian of basis elements in kerqg", does not have zeroes. 
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From the intertwining condition (|63p . we derive the following system of equations: 

V2{x)-Vi{x)^2a'{x), (67) 

a{x)[V2{x) - Viix)] - Wlix) = a"ix) + 2p'{x), (68) 

P{x)[V2{x) - Vi{x)] - 2aix)Vl{x) - Wl'{x) = /?"(a;) + 27'(a;), (69) 

-/ix)[V2{x) ~ Vi{x)] ~ (3ix)Vl{x) ~ a{x)Vl'{x) - Vl"{x) = Yix). (70) 

Relations (f66|) . (I67p and condition 1 imply that V2(a;) G C'g~^- We deduce from ((66|) . ([68|) . condition 1, and 
the inclusion V2(a;) G C^^^ that Q;"(a:;) G C'k"^! *-e-j Qf(x) G Cg. From ([67|) . condition 1 and the fact that 
Q;(a;) G C" it follows that ^2(2;) G C^^^ . It follows from ^E^j, (HH), condition 1 and the inclusions V2{x) G C^"^ 
and a{x) G that (3{x) G C^. From ((68|) . condition 1, and the inclusions V2{x) G C^~^, and Q;(a;) and 
l3{x) G it follows that a{x) G C'^''"^. Finally, it follows from ([57)) . condition 1 and the inclusion a{x) G Cj^'*'"'^ 
that V2(a;) G C^. Thus, the first statement is proved. 

The validity of equality is easily verified with the help of straightforward calculations. The fact that 
g2ix) and 50(2;) are in C^^^ and C^~^ , respectively, and equality (p5|) are obvious. Finally, to show that gi{x) 
belongs to C^'^^ we refer to the equalities 

3Viix) + a'{x) + 2(3{x) gC;1+\ (71) 

3V;{x)+p'{x) + 2j{x)eC^, (72) 

and 

Vlix)+j{x)eC^, (73) 

which, in turn, follow from equalities (l68l)-(l2Ql) since Vi,2{x) G C^, a{x) G C^"^\ /3(a;) G and 7(0;) G C^^^ 
Lemma 1 is proved. 

Corollary 3. By calculating the coefficient at in with the help of ([32|) . ((33|) and (|42|) . we deduce that 

a{x) = g2{x) = --^^ = -wi{x). (74) 

Hence, under the conditions of Lemma 1 

wi{x) G Cl+\ and VFi(a;) = Ce^''^^^")'^" G (75) 

8. Parametric formulas for coefficients of intertwining operators 

It was shown in [M] that the potentials Vi{x), V2{x) and coefficients of the intertwining operator can be 
parameterized by a single function which was denoted W{x) in [M]. It is not difficult to check that this function 
is connected with G{x) by the relation 

I^ = G-i(Ai+A2 + A3). (76) 
For convenience, we give parametric formulas obtained in [34 in the notation of the present work. 

Via = al T .92 - 2G + Ai + A2 + A3, (77) 

^1 



5i = :^[52-3G + Ai + A2 + A3], (78) 



1 



50 = al - .92 [.92 - 3G + Ai + A2 + A3] + - V(G')2+4P3(G). (79) 

Let us emphasize that, in contrast to ([5^ . parameterizations ([77|) - ([75|) are valid for any case. Parameterization 
([55)1 which supplements ([77)1 - ([75]) . (see, in addition, ([Ti]) ) is derived in [M] as well, but only under the following 
conditions 

G{x) + Const, \G'{x)f + 4P3(G(a;)) ^ 0, (80) 

whereas in the present work we derive this parameterization under a weaker condition (|58p . The case of 
conditions ([50]) , (|6ip was not considered in ^34j . 
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9. Relations between parameterization function and partial Wron- 
skians 

In the following lemma, we indicate basic relations between the parameterization function G{x) and Wronskians 
of a part of canonical basis elements in the kernel of . 

Lemma 2. Assume that: 1) the conditions of Lemma 1 are fulfilled; 2) Wjk — (p'jCf'k—'Pj'P'k^ 3) \/ (G'Y + '^PziG) 
is the same branch of the root as above. Then: 

1) if a Jordan form of the matrix S for the operator contains three Jordan cells of first order, 

h+(f>i = Ai0i, /l+02 = \2(f>2, h+(f)3 = A303, (81) 

then the following identities hold: 



2 ' 



(82) 



G - A, - £,,,(A, - A,)(A, - Xi)^^, (83) 



G' + V(G')2+4P3(G) = 2(Ai - A2)(A2 - A3)(A3 - Ai)^^^i^, (84) 

W12W23W31 



G' - V(G')^+4P3(G) - 2(Ai - A2)(A2 - A3)(A3 - Ai) '\ (85) 

where j,k,l is an arbitrary permutation of 1,2,3 and summation for a repeated index is not performed; in 
addition, the branch of \J (G')^ + 4P3(G) is independent of the choice of numbering of canonical basis elements; 

2) if a Jordan form of the matrix S for the operator contains a Jordan cell of first order and a Jordan 
cell of second order: 

= Ai0i, /l+02 = A202 + 03, ^^03=A303, A3 = A2, (86) 

then the following identities hold: 



Wi, = iX2-M)Ms, W^, = -<t>l = (^)' = (A 



A2)— (87) 



G-Ai = (Ai-A2)"^::^^, G-A2 = (Ai -A2)- 



a2 



G' + yW+4JMG) = -2(Ai - A2)^^, G' - VW+Tf^G) = 2(Ai - \2? ^^^^^\ (89) 



and the branch of ^ (G')^ + 4P3(G) is independent of the choice of numbering of canonical basis elements; 
3) if the Jordan cell of the matrix S for the operator q^ consists of a single Jordan cell of third order: 

h^(j)l = Xl(f>l + (t)2, h^(f>2 ^ X2(j)2 + (1)3, h^4'3 ^ X3<j)3, A3 = A2 = Ai, (90) 

then the following identities hold: 



(91) 



G-A,.^, (92) 

G' + V(G')^+4P3(G) = -2|^, G'- V(G')'+4P3(G) = 2^. (93) 

Wi W-^ 

Proof. Identities ((82l) . ((87l) and (jOT]) are easily checked with the help of straightforward calculations in 
which we use relations jM]) and ([90| . 
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Identities ([55]) . ([55)1 and ([5^ follow from identity (|54p - for convenience, we write the latter identity in the 
case considered in the form 

G-A. = -^(^^; (94) 



as well, those identities follow from identities ((82|) . ((87)) . ([91]) and from the fact that we can renumber elements 
of canonical basis so that any given eigenvalue of the matrix S for the operator gets index 2 (sec (|94|) ). 

In the case G ^ , j 1, 2, 3 identities (Ull), ((851), dM]), and ([Ml) follow from identities ([83 l) . (jSSl and ([92]) 
and relations dM]), ([SIl). Before we prove identities jMl), (|H5|), dM]) and dM]) in the case G = Aj, j = 1,2,3 let 
us show that this case is equivalent to the existence of two Jordan cells with the same eigenvalue Xj in a Jordan 
form of the matrix S for the operator . In the latter case, the validity of the desired identities is obvious. 

Identities d82|) . d87| . |9T|l and the fact, that any formal eigenfunction of the Hamiltonian (different from 
identical zero) can have zeroes of first order only imply that the right-hand sides of expressions d83| . d88| and 
d92]) either are identical zeroes (if there are two Jordan cells with the same eigenvalue in a Jordan form of the 
matrix S for the operator q'^) or can have zeroes of the order not exceeding four. Thus, the identity G = Aj 
holds on some interval if and only if this identity holds on whole axis which is equivalent to the existence of 
two Jordan cells with the same eigenvalue \j in a Jordan form of the matrix S for the operator q^ . Hence, 
identities dUl), dSSl), dSSl) and (l93|) are proved. 

To show that the branch of -y/ (G')^ + 4P3(G) is independent of the numbering of canonical basis elements, 
one can renumber these elements, derive formulas similar to (|84p . (|85p . (j89p and (|93p for new numbering and 
compare the results. Lemma 2 is proved. 

Corollary 4. In the proof of Lemma 2, it was shown that for any j the function G{x) — Xj either is identical 
zero on the whole axis or can have zeroes only of order not exceeding four. In addition, the relation G{x) = Xj 
is equivalent to the existence of two Jordan cells with the same eigenvalue Xj in a Jordan form of the matrix S 
for the operator q^ , which in view of Theorem 2, is equivalent to the possibility to strip-off the operators q^ . 

Corollary 5. Under the conditions of Lemma 1, the functions (/>i(a;), (f>2{x) and (j)3{x) belong to G^^^; 
hence it follows from identities dH]), d87)) and (|9T|) . that the Wronskians Wki{z) (case 1), Wi3{x) and W23(x) 
(case 2), and W2{x) (case 3) belong to C^'^^. We apply these inclusions, take differences of identities (l84|) and 
d85|) : |89|) : |93|) . refer to inclusions d75|) . and take into account that Wi{x) has no zeroes (as the Wronskian of 



a basis in ker q^ ) to show that 



V{G'y+4Ps{G)eC^+'. (95) 
In its turn, inclusion dSS]) and formulas dSi)) . dM]) and (1^ provide that 

G{x) G G^'. (96) 

Corollary 6. If coefficients of q^ are real, then coefficients of the polynomial P^i^h^) = qfqf are real as 
well. Hence, either all of the numbers Xj are real or one of these numbers is real and two are mutually complex 
conjugate. Without loss of generality, we assume that elements of the canonical basis in ker^g" that correspond 
to real A^-, are chosen real, and elements that correspond to complex conjugate Xj are complex conjugate. Then, 
if all of Xj are real, the root ^(G')^ + 4P3(G) is real for any a; e R in view of dMl), dM]) and dM]). If there is pair 
of complex conjugate values Xj (obviously, this is possible only if all of Xj are different), then Wi{x) is purely 
imaginary (since complex conjugation of Wi{x) corresponds to a permutation of two lines in the definition of 
Wi{x)) and ^/((F)^ + TP^(G) is also real for any a; e M by virtue of dUl). Thus, 

[G'(x)]2 + iP3{G{x)) ^0, xeR. (97) 



10. Lower bound of the parameterization function 

A lower bound for the parameterization function G{x) is given by the following lemma. 

Lemma 3. //, under the conditions of Lemma 1, the intertwining operator q^ cannot he stripped-off, 
coefficients of q^ are real and A3 is the minimal real eigenvalue of the matrix S for the operator q^ , then the 
inequality 

G{x) > A3, xe M, (98) 
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holds. Proof. First we show that the inequahty 



G{x) ^ A3, x£R (99) 

takes place. Assume that for some pomt xq E M, inequahty is violated. Then P3{G{xo)) < and, 
consequently, by (P7)) ,the derivative G'{xo) / 0. Moreover, G'{x) does not vanish on the entire interval 
which contains xq and on which the inequality G{x) < A3 holds. Hence, G{x) either strictly increases or 
strictly decreases on this interval. Obviously the interval is not bounded from the left (right), if G{x) increases 
(decreases) on it. Let us show that the assumption about the violation of (|99p leads to a contradiction. For 
definiteness, we consider the case where G{x) increases on the above-mentioned interval. The case of decreasing 
G{x) is treated similarly. By inequality ([W]) . the inequality G'{x)/y/—P3{G{x)) ^ 2 is valid for any point of 
the considered interval. Integrating the latter inequality from x to xq, we deduce that 

G{xo) 

dG 

—==^2{xo-x), x<xo. (100) 

The left-hand side of inequality pOOp is bounded for x —00 while its right-hand side tends to -f 00. This 
contradiction proofs inequality ([M)) . 

To prove that G{x) — A3 has no zeroes, we use identity ([55]) which expresses ^2(2;) in terms of G{x) (see 
also ([7^ V Let us assume that there is a point € K such that G{xo) = A3. Since cannot be stripped-off, 
Corollary 4 and inequality l|99p at the point xq imply that the function G{x) — A3 has a zero of even order 2n, 
G'{x) has a zero of order 2n — 1, and G" [x) has a zero of order 2n — 2, where n is either 1 or 2; in addition, it 
is obvious that 

G"{xo) ^ 0. (101) 

First we consider the case where A3 is a zero of -P3(A) of order one. In this case, inequality ([M)l and the 
fact that the order of the root of G{x) — A3 is even imply condition ((58)) . which allows us to use formula ((59)) . 
Finally, since 

P^(G(xo)) = (A3 - Ai)(A3 " A2) > 0, (102) 

and inequality (|10ip holds, the right-hand side of ([5^ at the point xq is infinite, which contradicts to (1751) . 
Hence, G{x) cannot equal A3. 

Now we assume that A3 is a zero of P3(A) of order two or three. In this case, the numerator of ((59l) has, 
obviously, a zero of order 2n — 2 at the point xq and the denominator has a zero of order 2n — 1. Hence, 52(2;) 
has a pole at the point cco, which is impossible. Thus, G{x) cannot equal A3, and Lemma 3 is proved. 



11. Theorem on reducibility of intertwining operators of the third 
order 

The assertion that any intertwining operator of the third order with real coefficients is really reducible is de- 
scribed by the following theorem. 

Theorem 5. Assume that the conditions of Lemma 1 are satisfied and that the intertwining operator cannot 
he stripped-off has real coefficients. Let A3 he the minimal real eigenvalue of the matrix S for the operator . 
Then there exist intertwining operators pf and kf^ of the first order andp^ and /c^ of the second order such that: 

1) coefficients of pf and are real, and, in addition, coefficients of these operators at helong to C^^^; 

2) coefficients of P2 and fc^ are real, and, in addition, coefficients of these operators at d and belong to 
Cg^^ and C^, respectively: 

3) 

pt^ipTY, k+^ik-y, pt^{P2)\ k+^ik^Y; (103) 
4 ) the matrices S for the operators pf and k^ consist of A3 
5) 

% ^k^Pi ^k^P2, qt "^^Pikt ^Ptkt, (104) 
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= hip^ , k2 hi = h k2 , P2 ~ ^2^2 ; k^h2 = h . (105) 
h+pi = p^hi, hik} = k^h^, h+p^ = p^h2, h2kt = k^h^, (106) 

where hi and /12 are intermediate Hamiltonians with real potentials from . 

Proof. We consider the case with and fc^ only since the statements of Theorem 5 for the case of pf 
and k^ are easily verifiable with the help of transposition, and the statement for the cases of kf and pf follows 
from the symmetry between h~^ and h~ . 

Let us define p^ and k^ by the equalities 

Pi^r^, k2=rir2. (107) 

Then existence of an intermediate Hamiltonian hi and intertwining (|105|) follows from relations ((35)) and p8|) . 
The fact that the potential of the Hamiltonian hi given by the formula 

Vi{x) - 2[ln03(.T)]" = Vi{x) - 2w'j,{x) (108) 

(see (HU) is real and belongs to the space follows from ((56|) . (f74|) . ([97]) . (fOS]) . from the fact that (72 (2^) is real, 
and from the inclusions ([M)) . and .92(2^) G C'r^^ (see Lemma 1). In addition, the function 

X3(x) = = -U.3, (109) 

which is the coefiicient at 9° of the operator pj~ , is obviously real and belongs to Cj^^^ . To prove that coefficients 
of k2 are real and belong to the spaces of smooth functions of Theorem 5 we first apply relations ((33| . (l42l) . 
([5T|l. ((52)) and ([54]) to transform k^ to the form 

A:2" =d^ + [w-i - wi)d +{G + Vi-wl~ W1W3 ~ 2A3), (110) 

and then take into account the following statements: Vi{x), wi{x) = —32(2;), A3 (see the Theorem 5 conditions), 
G{x) and Wsix) are real, identity ^ and inclusions dM]), and ^ hold, Wsix) € C2+\ and Vi{x) e C^. 
Finally, the fact that the matrix S for the operator p^ consists of A3 follows from psp . Theorem 5 is proved. 
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